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ABSTRACT 

This  paper  deals  with  the  problem  of  determining  the  economic  order  quantity  (EOQ),  as  a  function  of  the  setup 
cost  and  the  holding  cost  in  the  interval  sense.  Practically  vagueness  caused  by  the  variation  in  fixing  these  costs  is 
inevitable.  Intuitionistic  fuzzy  inventory  model  with  instantaneous  replenishment  and  no  shortages  is  analyzed  to  compute 
the  economic  order  quantity  and  the  total  annual  cost  by  assigning  fuzzy  quantity  and  intuitionistic  fuzzy  quantity  instead 
of  real  quantity  to  these  costs.  Parametric  programming  technique  is  applied  and  the  results  are  compared  numerically  both 
in  fuzzy  optimization  and  intuitionistic  fuzzy  optimization  techniques.  Necessary  graphical  presentations  are  also  given 
besides  numerical  illustrations. 
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1  INTRODUCTION 

Inventory  problems  are  common  in  manufacturing,  maintenance  cost  and  business  activities  in  general.  Recently 
much  attention  has  been  focused  on  EOQ  models  with  fuzzy  carrying  cost,  fuzzy  shortage  cost,  fuzzy  setup  cost,  fuzzy 
demand  etc;  this  means  that  elements  of  carrying  cost,  shortage  cost,  setup  costs  and  demand  are  fuzzy  members 
[10], [20], [23].  Zimmerman  [28,29,30]  showed  the  classical  algorithms  can  be  used  in  few  inventory  models. 

Inventory  model  such  as  instantaneous  stock  replenishment  and  no  shortages  is  analyzed  by  assigning  fuzzy 
quantities  to  the  setup  costs  (Cs)  and  holding  cost  (Q)  instead  of  crisp  values. 

Recently  fuzzy  concept  is  introduced  in  the  inventory  problems  by  several  researchers.  Park  [21],  Vujosevic[27] 
et.al,  Chang[5,7]  et.al,  Liu[19]  et.al  are  proposed  the  EOQ  model  in  the  fuzzy  sense  where  inventory  parameters  are 
triangular  fuzzy  concept  in  decision. 

An  early  work  using  fuzzy  concept  in  decision  making  has  been  performed  by  Bellman  and  Zadeh[4],  through 
introducing  fuzzy  goals,  costs  and  constraints.  Lee  et.al[18]  introduced  the  application  of  fuzzy  set  theory  to  lot  sizing  in 
material  requirements  planning.  In  their  paper  uncertainty  in  demand  is  modeled  by  using  triangular  fuzzy  numbers. 

One  of  the  interesting  generalization  of  the  theory  of  fuzzy  sets  is  the  theory  of  intuitionistic  fuzzy  sets  introduced 
by  Atanassov[l,16,17]  and  seems  to  be  applicable  to  real  world  problems.  The  concept  of  intuitionistic  fuzzy  sets  can  be 
viewed  as  an  available  information  is  not  sufficient  for  the  definition  of  an  imprecise.  Intuitionistic  fuzzy  sets  are  fuzzy  sets 
described  by  two  functions;  a  membership  function  and  a  non-membership  function  that  are  loosely  related. 
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Intuitionistic  fuzzy  set  can  be  used  to  simulate  human  decision-making  process  and  any  activities  requiring  human 
expertise  and  knowledge  that  are  inevitably  imprecise  or  totally  reliable.  An  interval  valued  intuitionistic  fuzzy  sets  are 
analyzed  by  Atanassov  and  Gargov[l  1,12,13]  Atanassov  and  Rreinorich[16,17]  implemented  intuitionistic  fuzzy 
interpretation  of  interval  data.  The  crisp  values  are  compared  numerically  both  in  fuzzy  optimization  and  intuitionistic 
fuzzy  optimization  techniques.  Objective  of  this  paper  is  to  find  Q°  and  TC°  in  both  intuitionistic  fuzzy  optimization 
method  and  fuzzy  optimization  method.  Necessary  graphical  presentations  and  numerical  illustrations  are  also  given. 

2  FORMULATION  OF  PARAMETRIC  PROGRAMMING  PROBLEM  FOR  SOME  BASIC 
INVENTORY  MODELS 

The  setup  cost  X,  holding  cost  Y±,  and  shortage  cost  Y2  are  approximately  known  and  are  represented  by  the 
following  intuitionistic  fuzzy  sets: 

X  =  {Qc,nx(x),vg(x))\xeX)  (2.1) 
%  =  {(yum  (yd,  v?i  (yd)  I  Vi  e  YJ,  (i  =  1,2)  (2.2) 

where  X,  Y±  and  Y2  are  crisp  universal  sets  of  setup  cost,  holding  cost  and  shortage  cost  respectively  and  fJ-x(x) 
and  Myj(yj),  (t  =  1,2)  are  the  respective  membership  functions. 

Let  a  denotes  the  minimal  acceptable  degree  of  objectives  and  constraints  and  /?  denotes  the  maximal  degree  of 
rejection  of  objectives  and  constraints.  Then  the  (a,  /?)-level  of  X,  Yv  (i  —  1,2),  [6,10,1 1],  are 


Xa  =  {(x  EX\fix>  a)},  (2.3) 

Yia  =  [(yi£Yi\(in>cc)},(i  =  1,2)  (2.4) 

XP  =  {(x  e  X\vx  <  /?)}  (2.5) 

Yi"  =  {(yt  6  Yt\v?i  <p)},  (i  =  1,2)  (2.6) 


where  a  >  ft;  a  +  /?  <  1  and  a,  ft  >  0. 

The  quantities  Xa,  Yia  (i  —  1,2),  X^  and  Y^,(i  —  1,2)  are  crisp  sets.  Using  (a,/?)-level  the  setup  cost,  holding 
cost  and  shortage  cost  can  be  represented  by  different  levels  of  confidence  intervals  [16,17,25,26].  Hence  an  intuitionistic 
fuzzy  inventory  models  can  be  reduced  to  a  family  of  crisp  inventory  models,  with  different  {a,  /?)-level  cuts 

{Xa\0<a<  1},  {Yla\  0  <  a  <  1},  (i  =  1,2), 
{Xp\0  <  p  <  1}  and  {Ytp\Q  <  p  <  l}(i  =  1,2) 

The  above  sets  represent  sets  of  movable  boundaries  and  they  form  nested  structure  for  expressing  the  relationship 
between  the  crisp  sets  and  intuitionistic  fuzzy  sets. 

Let  the  confidence  intervals  of  the  intuitionistic  fuzzy  sets  X  and  Yh  (i  =  1,2)  be  [lXcc>uxa]  and  Uyt  >uyi  L 
(i  —  1,2),  [uxp,lxp]  and  [uy.p,  ly.p],  (i  —  1,2)  respectively.  Since  the  set-up  cost,  holding  cost  and  shortage  cost  are 
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intuitionistic  fuzzy  numbers,  using  Atanassov's  extension  principle  [16,17,25],  the  membership  and  non-membership 
functions  of  the  performance  measure  p(X  and  5^),  (i  =  1,2)  are  defined  as 

MP(*,n)(z)  =  suP*ex,y;en  miri  0^00^60/*  =  P(x,yd},(i  =  1,2)  (2.7) 
and 

VpixfjOO  =  infxex,yien  max  {vg(x), v^Cy^/z  =  p(x,yj},(i  =  1,2)  (2.8) 

Construction  of  the  membership  function  Mp(x,fi)(z)'   (i  =  1,2)  and  vp(^  j?j(z)(i  =  1,2),  are  respectively 
equivalent  to  say  that  derivation  of  (a,  (3)  -  levels  of  /Wj^f-jCz),  (i  =  1,2)  and  vp(jj>jf)(z) 

(i  =  1,2), 

From  the  equation  (2.7)  the  equations  /ip(^  j?)(z)  =  a,  (i  =  1,2)  is  true  only  when  Hx(x)  —  a, 
^(yd  ^  a  or  WW  ^  «»  ^Oi)  =  «  are  true. 

From  the  equation  (2.8)  the  equations  vp(^f)(z)  =  /?,  (i  =  1,2)  are  true  only  when  Vjj(;t)  =  /?, 

v^(yf)  <  /?  or  Vxix)  <  /?,  v7l{yd  =  P  are  true. 

The  parametric  programming  problems  have  the  following  form: 

lVa  =  minp(x,yi)  (2.9) 
such  that 
L  <x<ur, 

lyia  ^Vi^  uyia>  0  =  W). 
and 

uPa  =  max  p(x,y0 
such  that 

lXa<x<  uXa,  (2.10) 
'yia  ^  yt  ^  "yta.  G  =  1.2), 

upp  =  min  p(x,y;)  (2.11) 
such  that 

uxp  <x  <  lxp, 

uy,fs  <yt<  ly.fs,  (i  =  1,2), 
and 

lpp  =  maxp(x,yi) 
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such  that 

uxp  <  x  <  lxn,  (2.12) 
uy.n  <yt<  ly.ft,  (i  =  1,2), 

If  both  lPa  and  uVa  are  invertible  with  respect  to  a  then  the  left  shape  function 

L(z)  =  l~1pa  and  the  right  shape  function  R(z)  —  u~1pa  [2,3,24]  can  be  obtained.  From  this  the  membership 
function  fJ-p^x.Y^i2)^  0  =  1<2)  is  constructed  as 

!L(z)  for  z1  <  z  <  z2 
\forz2<z<z3  (2.13) 
R  (z)  for  z3  <  z  <  z4 

where  zr  <  z2  <  z3  <  z4,  L^)  =  R(z4)  =  0  and  L(z2)  =  R(z3)  =  1. 

If  both  up/?  and  lpp  are  invertible  with  respect  to  /?  then  the  left  shape  function 

L(z)  =  u_1p^  and  the  right  shape  function  R{z)  —  /_1p^  [2,3,24]  can  be  obtained.  From  which  the  membership 
function  vp^  j?)(z),  (t  =  1,2)  is  constructed  as 

!1  —  L(z)  for  z1<  z  <  z2 
0forz2<z<z3  (2.14) 
1  —  R(z)  for  z3  <  z  <  z4 

where  z-l<z2<z3<  z4,  LCzj^)  =  fi(z4)  =  1  and  L(z2)  =  R(z3)  =  0. 

3.  MODEL  -  EOQ  PROBLEMS  WITH  INSTANTANEOUS  REPLENISHMENT  AND  NO  SHORTAGES 
Assumptions 

3.1.  The  inventory  system  pertains  to  a  single  item. 

3.2  Annual  Demand  (D)  is  deterministic. 

3.3  The  inventory  is  replenished  in  a  single  delivery  for  each  order. 

3.4  Replenishment  is  instantaneous. 

3.5  There  is  no  lead  time. 

3.6  Shortages  are  not  allowed. 

Using  the  concept  of  a,  ft  level,  the  above  intuitionistic  fuzzy  inventory  models  can  be  reduced  as  EOQ  model 
with  instantaneous  replenishment  and  no  shortage  [9,14,15,22]  for  which 

<?°  =  Pgf  (3.1) 
and 

TC°  =f±]  +  [^]  =  [2DCsC1]l  (3.2) 
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where  Cs  and  C1  represent  the  set-up  cost  and  the  holding  cost  respectively. 
4.  ILLUSTRATIONS 

Consider  an  EOQ  problem  with  instantaneous  replenishment,  without  shortages.  The  setup  cost  and  holding  cost 
are  fuzzy  numbers  represented  by  X  =[400,450,550,600]  and  J^=[0.7 ,0.8,1. 1,1.1].  The  a  level  of  the  membership 
functions  JUe(x)  ,  jU^(yi))sre  [400+50a,600-50a]  and  [0.7+0. la,  1.1 -0.1  a]  respectively.  The  (3  level  of  the 
non  membership  functions  V^(x)  ,  VTi  (yi)  are  [450-500,  550+50/?]  and  [0.8-0.1  0,1+0.1  /?]  respectively. 

From  the  equation  (2.9),  (2.10),  (2.11)  and  (2.12),  the  parametric  programming  problems  are  formulated  to  derive 
the  membership  function  for  Q 

They  are  of  the  form 


(3.3) 


with  400+  50a  <x< 600- 50a 


0.7+0.1  a  <y,  <  1.1  -0.1a 


(3.4) 


with  400+  50a  <x< 600- 50a 


0.7+0.1  a  <y,  <  1.1  -0.1a 


Where  0<a<  1 


(3.5) 


with  450  -  50/?  <  x  <  550  +  50/5 


0.8-0.1  P  <yi  <1  +  0.1/7 


(3.6) 


with  450  -  50/?  <  x  <  550  +  50/? 


0.8-0.1  /?  <y,  <1  +  0.1/3 


Where  0  <  /?  <  1 
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^0  is  found  when  x  and  y{  approach  their  lower  and  higher  bound  respectively.  Taking  D=  10,000  units, 
the  optimal  solution  for  (3.3)  is 


lQl 


20xl04x(400  +  50a) 


-il/2 


(3.7) 


Also  Wg0  is  found  when  x  and  yj  approach  their  upper  bound  and  lower  bound  respectively.  In  this  case  the 
optimal  solution  for  (3.4)  is 


UQl  = 


20xl04x(600  -  50a) 
1  +  a 


1/2 


(3.8) 


/  o/}  is  found  when  x  and  yj  approach  their  higher  and  lower  bound  respectively.  Taking  D=  10,000  units,  the 
optimal  solution  for  (3.5)  is 


I   P  = 

a* 


20xl04x(450  -  50/?) 
10  +  /? 


1/2 


(3.9) 


Also  u  0iS  is  found  when  x  and  yi  approach  their  lower  bound  and  higher  bound  respectively.  In  this  case  the 
optimal  solution  for  (3.6)  is 


20xl04x(550  +  50/?) 
8-/? 


1/2 


(3.10) 


The  membership  function  jU—(z)  is  obtained  and  given  by 


r~Qo 


1 1  z2  —  Rx  1 07 

,       7     for  2696.7995 <  z  <  3000 
z2  +  107 

1  for  3000  <  z  <  3708.0992 

12xl07  —  772 

 " — P-  for  3708.0992  <  z  <  4140.3934 


10'  +zl 


(3.11) 


The  graph  of  JU—(z)  is  depicted  in  Figure  1 . 


The  membership  function  V— (z)  is  obtained  and  given  by 
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—2 — for  2696.7995<  z  <  3000 


z'+lO' 


for  3000<z<  3708.0992 


(3.12) 


8z2-llxl07 


z2  +  l0 


Y —  for  3708.0992  <  z  <  4140.3934 


The  graph  of  ^^(-Z)  is  depicted  in  Figure  2. 

The  parametric  programming  problem  corresponding  to  the  total  annual  cost  TC°  has  different  only  from  the 
(3.3)  and  (3.4)  in  the  objective  function  and  given  below 


/„„„  =min{2DxyiJ 

'I  a 


(3.13) 


and  U   0    =max  {2Dxyi] 

T C  a 


(3.14) 


The  parametric  programming  problem  corresponding  to  the  total  annual  cost  TC°has  different  only  from  the  (3.5) 
and  (3.6)  in  the  objective  function  and  given  below 


I    g  =max  {2Dxy!}1/2 

TC° 


and  U     a  =min  {2Dxy!j 

TC" 


(3.15) 
(3.16) 


From  the  above  problems  (3.13)  &  (3.14)  I    ,  is  obtained  when  both  x  and  yi  approach  their  lower  bound  U  „ 

TC   a  TC  a 

is  obtained  when  both  x  and  yi  approach  their  upper  bound.  Taking  D=10,000  units. 


/     =[2xl04x(280  +  75a  +  5a/)]: 


(3.17) 


2\12 


u  0  =[2xl04x(660-115a  +  5az)] 

TC  a 


(3.18) 


The  membership  function  jU    „  (z)  is  given  below 


-150xl02+Vl06  +  40z2 
2x1b1 

I 

23xl03-Vl06+40z2 


for  2366.4319  <  z  <  2683.2816 
for  2683.2816<z<3316.6248 


2xlOj 


for  3316.6248  <  z  <  3633.1804 


(3.19) 
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The  graph  of  jU—(z)  is  given  in  Figure  3.  For  the  crisp  values  Cs=  ?  500, Cl=    ?  0.9,  D=10,000  units 
Q°=  3333.333  units  and  TC°=  ?  3000. 

But  defuzzification  [6,8]  of  jU—(z)  and  jU—(z)  give  Q°=  3354.05  units  and  TC°  =  ?  2999.9532  respectively. 


So  if  definite  fluctuation  of    ?  50/-  and  the  maximum  possible  fluctuation  of    ?  100/-  in  Csand  definite 

fluctuation  of  ?  0.1  and  the  maximum  possible  fluctuation  of  ?  0.2  in  Cl  are  identified,  then  the  amount  of  Q°  to  be 
increase  suitably  (it  can  be  computed  by  using  the  above  techniques)  so  as  to  maintain  the  total  annual  cost. 

The  graph  of  V—(z)  is  given  in  Figure  4. 


But  defuzzification  [6,8]  of  V—(z)  and  V—(z)  give  Q°=  3354.05  units  and  TC°  =  t  2999.9532  respectively. 

So  if  definite  fluctuation  of    T  50/-  and  the  maximum  possible  fluctuation  of    ?  100/-  in  Csand  definite 

fluctuation  of  ?  0.1  and  the  maximum  possible  fluctuation  of  ?  0.2  in  Cl  are  identified,  then  the  amount  of  Q°  to  be 
increase  suitably  (it  can  be  computed  by  using  the  above  technique)  so  as  to  maintain  the  total  annual  cost. 

From  the  above  problem  (3.15)  &  (3.16)  /     »  is  obtained  when  both  x  and  yl  approach  their  upper  bound 

TC" 

U     a  is  obtained  when  both  x  and  y[  approach  their  lower  bond.  Taking  D=  10,000  units, 

TC" 

1 

I    ,  =  [2xl04  x  (360 -S5J3  + 5  J32)]2  (3.20) 

TC" 


u    ,  =[2xl04x(550  +  105/?  +  5/?2)]2 

TC" 

(3.21) 


The  non  membership  function  V —     is  given  below 

TC"(Z) 


-3  for  2366.4319  <z<  2683.2816 


v — 

TC"(Z) 


CONCLUSIONS 


2x10' 

1  for  2683.2816<z<3316.6248  (3.22) 

-21xl03  +  Vl06  +  40z2 


2xlOj 


for  3316.6248  <  z  <  3633.1804 


In  this  paper,  the  membership  functions  and  the  non  membership  functions  for  the  economic  order  quantity 
Q"  and  the  total  cost  TC"  axe  obtained  for  the  elementary  intuitionistic  fuzzy  inventory  models. 
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